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EQUILIBRIUM SPECTRA OF SECONDARY COSMIC-RAY 


POSITRONS IN THE GALAXY AND THE SPECTRUM 
OF COSMIC- GAMMA-RAYS RESULTING FROM THEIR ANNIHILATION 

F. W. Stecker 
ABSTRACT 

A general formula is derived for calculating the gamma-ray spectrum re- 
sulting from the annihilation of cosmic-ray positrons. This formula is used to 
calculate annihilation-gamma-ray spectra from various equilibrium spectra of 
secondary galactic positrons. These spectra are then compared with the gamma- 
ray spectra produced by other astrophysical processes. 

Particular attention is payed to the form of the gamma-ray spectrum re- 
sulting from the annihilation of positrons having kinetic energies below 5 keV. 
It is found that for mean leakage times out of the galaxy of less than 400 million 
years, most of the positrons annihilating near rest come from the beta decay of 
unstable nuclei produced in cosmic-ray p-C 12 ,p-N 14 , and p-O 16 interactions, 
rather than from pi-meson decay. It is further found that the large majority of 
these positrons will annihilate from an S state of positronium and that 3/4 of 
these will produce a three-photon annihilation continuum rather than the two- 
photon line spectrum at 0.51 MeV. The results of numerical calculations of the 
gamma-ray fluxes from these processes are given. 
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EQUILIBRIUM SPECTRA OF SECONDARY COSMIC-RAY 


POSITRONS IN THE GALAXY AND THE SPECTRUM 
OF COSMIC- GAMMA-RAYS RESULTING FROM THEIR ANNIHILATION 

F. W. Stecker* 

Goddard Space Flight Center 
Greenbelt, Maryland 

I. INTRODUCTION 

The annihilation of cosmic-ray positrons has for some time been recognized 
as a potential source of cosmic gamma rays. Gamma-ray fluxes from cosmic 
positron annihilation have been estimated and discussed by various authors 
(Pollack and Fazio, 1963; Hayakawa, et. al., 1964; Ginzburg and Syrovatskii, 
1964a, b). Pollack and Fazio have discussed the possible relationship between 
the present flux of 0.5 MeV gamma-radiation from positron annihilation and the 
cosmic-ray intensity and galactic gas density 10 9 years ago. Ginzburg and 
Syrovatskii have pointed out that the intensity of the 0.5 MeV line may be a 
sensitive measure of the leakage rate of cosmic-ray positrons from the galaxy. 
These authors have also given an approximate formula for the calculation of the 
gamma- ray spectrum from the annihilation of high-energy positrons. 

It has therefore become apparent that the cosmic-gamma ray spectrum 
from cosmic- ray positron annihilation may contain much potential information 
reflecting various astrophysical conditions in our own galaxy and in possible 

•Research Associate— National Academy of Sciences— National Research Council. 
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cosmic-ray sources, both galactic and extragalactic . It is for this reason, that 
the author had recently undertaken a more detailed investigation of the various 
aspects of the cosmic positron annihilation problem (Stecker, 1967b). That work 
has now lead to this further treatment pointing out the potential richness and 
complexity of the problem. 

The point of departure for our discussion will be a derivation of the general 
formula for calculating the annihilation- gamma-ray spectrum (AGS) from cosmic- 
ray positrons valid at all energies. This formula is then used to calculate AGS 
from various equilibrium spectra of secondary cosmic- ray positrons. The 
results of various numerical calculations will be given and discussed. 

The significance of direct measurement of the positron equilibrium spectrum 
itself will be discussed as a possible sensitive indicator of the diffusion and 
leakage of cosmic-ray leptons out of the galaxy. 

II. THE ANNIHILATION- GAMMA- RAY SPECTRUM 

The cross section for positron annihilation as a function of energy was first 
deduced by Dirac (1930). An excellent presentation of the theory is given by Heitler 
(1960). The most important annihilation mode of the free electron-positron system 
is the annihilation 

e + + e~ -* y + y (1) 

The frequency of this annihilation is 372 times greater than that of the free 
three -photon annihilation and we may neglect all but the two-photon mode in 
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considering the gamma- ray spectrum from free e + - e“ annihilation. (However, 
as we shall see later, the three-photon mode becomes important when we consider 
the effect of positronium formation by positrons of energies less than 5 keV.) 

The differential cross section for gamma-ray production in the collision 
c.m.s. of a free two-photon annihilation may be written as 

, -o 

da - <Kxl y)dy (2) 

where cr 0 = ttt 0 2 , r 0 = e 2 /Mc 2 is the classical radius of the electron, y is the 
Larentz factor of the positron, /3 n Vl - l/y 2 , the c.m.s. Lorentz factor and 
velocity are given by 


r c - 311(1 fi e ■ /yrr ’ ( 3 ) 

X is the cosine of the angle between the incoming positron and the outgoing 
gamma- ray in the c.m.s., and the angular distribution function, <fi(x'> 7) is de- 
fined as 


<Kx\ 7) 


2/3*(i-yy 

i-eyy ' (i-/s c V ) 2 


(4) 


The energy of an annihilation-gamma- ray in the laboratory system is given by the 
Doppler relation as 


E y - (Me 2 ) y 2 (l + /3 C x) 


(5) 
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If we now define the dimensionless energy 77 = E y / Me 2 , we may use Eqs. 
(3) and (5) to determine the normalized distribution function uniquely relating 
X to the laboratory Lorentz factor of the positron and the laboratory energy of 
the gamma- ray as follows: 


f (x; 7. p) - 8 [x"x 0 ( 7?) ?>] 


where 


_ (27] -1) - y 

xo(v-r) - - ' jp . ~ 


( 8 ) 


The total number of gamma-rays produced per second by a positron of 
energy 7 Me 2 is 


Qtot ( 7 ) 


2 ff o n e^ c r 1 


27c 2 /3 C 


I 


dx 0 (y; 7) 


2°' 0 n e c r 1 


7 


£ 


dx<£(x; 7) 


where n is the number density of electrons in the medium. 
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Therefore, the gamma-ray source spectrum from the annihilation of cosmic- 
ray positrons with a density n(y) cm'V 1 is given by 


Q(v) 



(7) 


We now wish to integrate over the delta function in Eq. (7). The result of 
this integration is to replace the last integral in Eq. (7) by a function $( 77 , y) 
such that 

(i(xo 07. 7). y) for I xr 0 1 — 1 
7) = 1 . . ( 8 ) 

L 0 for W >J 

We may, therefore, replace the integral over dy by the algebraic function 
^(xo C 7 ?’ 7) > 7) provided the limiting condition on y 0 is transformed into a 
limiting condition on the integration over dy. The limiting conditions y 0 = ±1 
correspond to the relations 

V ± = 7 c 2 (l±£ c ) (9) 

It follows from Eq. (9) that the product 

’i.’i- = y (i-y) = y = ■ 2 4 i oo> 
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and the sum 


V + + V- ~ 27 2 -y + l 

Therefore, 77 + and 77. are the roots of the quadratic equation 

V 2 ~ 27 2 77 + 7 C 2 = 0 


which can be solved for y c in terms of 77, yielding 


y, 


2 

c 


77 2 7 + 1 

277-1 2 


(ii) 


( 12 ) 


(13) 


or 


_ 77 2 + (77 - l ) 2 

y 277-1 


(14) 


which is finite and positive for 77 >1/2 and has a minimum of y- 1 at 77 = 1. 

Figure 1 shows the curve defined by Eq. (14) split into the branch 
defined by the physical extremes y = ±1. Also shown are the asymptotes 
77 = 1/2 and 77 = 7 + 1/2, the line defining x = 0 and the shaded region corre- 
sponding to the physical values of I xl 1 1. It can be seen from Fig. 1 that the 
physical range of 7 defined by the shaded region (I xl 1 1) is bounded on the 
bottom by the curve of Eq. (14) and is unbounded on the top. 

Figure 1 indicates that no gamma rays can be produced from the annihila- 
tion process having energies less than or equal to Mc 2 /2. This physical 


6 


restriction may be seen more clearly as a direct consequence of Eq. (9) by 
noting that 


V- < V < v+ 


( 15 ) 


where 


1 - fir 


7 - = y 2 (l ~fi c ) - = TTT C > 1 


(16) 


so that ?)_ may approach, but never reach 1/2 as 7 - od. On the other hand 


v + - y 2 (l +£ c ) - 


(17) 


which increases without bound as y c 

The general restriction on the range of 17 may be designated by the intro- 
duction of the Heavyside step function, H + (rj 0 ] which is defined by the relation 


jl for V > Vo 

lo for 77 < r] 0 


We may therefore rewrite Eq. (7) in the form 


Q(v) 


4H + 



a 0 c 


(18) 


x 


f°° n(7 ) 

zjzrri $(xo( v,y),y) 

•'G(v) 1 


(19) 


7 


where 


G(t?) 


7? 2 + ( 7 ? ~ l ) 2 

2v - 1 


From the general results which we have obtained, we may immediately 
arrive at formulas for asymptotic spectra which may be used as guides in examin- 
ing the results of numerical calculations. These are obtained as follows: 

A) For two-photon annihilations at rest, it follows from Eq. (5) that the 
AGS is simply a line at energy 77 = 1 (0.51 MeV). This is, of course a familiar 
and expected result. 

B) The AGS from the two-photon annihilation of an ultrarelativistic posi- 
tron may be obtained from a consideration of the angular distribution function, 

<£(x; y) given by Eq. (4). At ultrarelativistic energies, the angular distribution 
function peaks sharply at y = ±1 so that the gamma- rays resulting from the 
annihilation lie close to the asymptotes of Fig. 1, viz. 17 = 1/2 and 77 = y + 1/2. 
This result, as pointed out by Heitler (op. cit.) may be discussed physically as 
follows. In an ultrarelativistic annihilation, the resulting photons are emitted 
in a sharply backward and sharply forward direction in the collision c.m.s. 
respectively. In the laboratory system, the forward photon carries off almost 
all the energy of the collision while the backward photon carries off an energy 
between about 0.25 and 0.5 MeV. Therefore, the AGS for 77 »1 may be obtained 
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by using the approximate production- cross- section 


( v < 7 ) - (7)8(17 -y) for 77 » 1 (20) 

The two-photon annihilation cross-section at ultrarelativistic energies has 
the asymptotic form 

O' 0 

(y) ~ ~[r& n (2y)-i] > y» 1 (2i) 

Therefore, the ultrarelativistic asymptotic form of the AGS is given by the 
expression 


n(7?) \ln{ 2 rj) - l] 

Q(V) ~ n e .(T 0 c ~ (22) 

as has been previously noted by Ginzburg and Syrovatskii. These asymptotic 
formulas should be kept in mind when examining the numerical evaluations of 
the AGS obtained from the exact formula given by Eq. (19). 

m. THE EQUILIBRIUM SPECTRUM OF SECONDARY GALACTIC POSITRONS 
FROM PI-MESON DECAY 

In order to calculate plausible annihilation-gamma- ray spectra from our own 
galaxy, it will be assumed that the only source of galactic cosmic-ray positrons 
above a few MeV energy is the result of primary cosmic-rays colliding with atoms 
of interstellar hydrogen and helium gas in the galaxy. These high-energy 
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collisions are known to produce positive pi-mesons which rapidly decay into 
positrons and neutrinos. Much accelerator data is available on the production 
of pi-mesons in interactions up to 30 GeV/c and many cosmic-ray studies have 
been made of higher energy interactions. Various calculations have been made 
using this data to calculate positron source spectra in the galaxy (Pollack and 
Fazio (1963); Ginzburg and Syrovatskii (1964b); Hayakawa, et al. (1964); Jones 
(1963); Ramaty and Lingenfelter (1966)). In these calculations, it is usually 
assumed that the primary cosmic- ray spectrum is essentially the same as that 
observed at the earth. Ramaty and Lingenfelter have shown that in the energy 
range where most of the pions are produced (above 500 MeV) the primary galactic 
cosmic-ray spectrum observed at the earth is little affected by local conditions 
in the solar system and that local effects induce very little error in the positron 
flux calculation. We therefore take for our source spectrum of cosmic-ray 
positrons, the one derived by Ramaty and Lingenfelter. (This spectrum is shown 
in Fig. 2.) 

We further assume that these positrons have reached a quasi-equilibrium 
condition in the galaxy determined by leakage out of the halo, annihilation, and 
energy loss by ionization and coulomb interactions, bremsstrahlung, Compton 
collisions and synchrotron radiation. The equilibrium density spectrum, n(y) , 
which will be used in Eq. (19), is then taken to be determined by the continuity 


equation 


where q + (y) is the positron source spectrum from pi-meson decay, r(y) is the 
total energy loss rate per unit time andr s (y) is the effective positron survival 
time given by 


1 _ 1 1 

t s (y) T \(y) + r i 00 


(24) 


with t a (y) and (y) being the annihilation time and mean leakage time for 
positrons of energy yMc 2 . These quantities are given by 


n 


(7) 


• °~o c 

y 


/F 


-^2 


y 2 + 4y + 1 


7 + 1 


L y 


2 _ 






(25) 


and 


1 _ /y 2 " 1 „ 

= yTT' T -e constant (26) 

where we have assumed a mean leakage time inversely proportional to velocity. 

Most of the interstellar gas in our galaxy is unionized with the exception of 
the so-called HII regions near the very hot O and B stars which are powerful 
sources of ionizing ultraviolet radiation. Allen (1963) gives the proportion of 
space near the galactic plane occupied by clouds of interstellar gas and dust as 
7% and that occupied by ionized clouds (HII regions) as 0.4%. We will therefore 
assume that the galactic gas is entirely neutral for the purposes of these 
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calculations and take for the energy loss rate from ionization the expression 


8 y 

r ! O') - T a o n e- 


22 + 4n {y(y- 1) (y 2 " l)} " 1.695 (^T~j ~ 


y 2 - 1 \ 1.39 

y 


(27) 


(See Heitler (1960); Morrison (1961)). 

The energy loss rate from bremsstrahlung may be taken as 


r B (y) = 7.3 x 10" 16 n e y (28) 

based on radiation lengths for hydrogen and helium given by Dovzhenko and 
Pomanskii (1964). 

The loss rate from synchrotron radiation and Compton collisions may be 
taken as 


r s+c (y) = 1.3 X 10' 9 (H 2 +3x 10 -11 p y )j 2 (29) 

where His given in gauss and the radiation density, p y , is given in ev/cm 3 
(Ramaty and Lingenfelter (1964)). 

The total positron energy-loss rate is taken to be the sum of Eqs. (27) - 
(29) and is simply denoted by r (y) . Equation (23) may then be solved to yield 
the equilibrium positron density spectrum in the form 

I CO 

dz q + (z) exp 


r 

Jry 


dw 


r(w)| t (w) 


( 30 ) 


(Stecker, (1967b)). 
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Equation (30) will be used to obtain numerical solutions for n(y), the posi- 
tron equilibrium flux 


i + (r) 


c n(-y) /y 2 - 1 

477 y 


(31) 


and the resulting AGS of 


Qiv) 


4H+ 



n cr n c 

e ~ 0 


4* 00 

*'G(7J) 


dry 


^(xp (jh r). r) 

yiy 1 ~ l I r(y)| 


f 00 

/* z 

X dz q y ( z ) exp 

Jy 

j dw 

J y |r(w)|r s (w) 


(32) 


In order to utilize Eqs. (30) - (32) to determine the positron equilibrium 
flux and AGS in the galaxy, we chose typical average values for the quantities 
n e > H, p y and Ty for both the galactic disc and the galactic halo. These values 
are given in Table 1. The radiation density, p , includes the contribution of 
0.25 eV/cm 3 from the 2.7°K universal microwave field (Stokes, Partridge and 
Wilkinson, 1967) 

Using the values given in Table 1 for the galactic halo to numerically in- 
tegrate Eqs. (30) - (32), typical galactic spectra are obtained for l + (y) and 
Q(p). The results are given in Figs. (3) and (4) for various values of mean 
path length X(g/cm 2 ) and corresponding Tf> The flux of annihilation-gamma-radiation 


13 



observed at the earth, I A ( 77 ) , is given by 


Ia(^> 



Q(v) 


(33) 


where L e f f is the effective path length for gamma-ray production. 

For leakage times less than 10 million years, the equilibrium positron 
flux has roughly the same characteristics of the source spectrum of Fig. 2 and 
its magnitude is proportional to the leakage time. In the case of longer leakage 
times, the positrons are trapped in the galaxy for a sufficient time for the energy 
loss processes, particularly ionization loss, to affect ;he spectrum by progressively 
flattening it below 30 MeV. Studies of spallating cosmic-ray nuclei yield a mean 
path length, X = 4 ± 1 g/cm 2 for cosmic-ray nuclei, corresponding to a mean 
leakage time of the order of 100 million years. The curve in Fig. 3 correspond- 
ing to X = 5 g/cm 2 is in agreement with the cosmic-ray positron measurements 
of Hartman (1967). It may be noted that measurements of the galactic positron 
flux below 30 MeV would yield a more sensitive check on the galactic mean leakage 
time. 

Figure 4 shows the annihilation-gamma-ray spectra obtained using the 
positron equilibrium fluxes of Fig. 3. The spectra shown are from annihilations- 
in-flight of positrons having energies greater than 5 keV. AGS from positrons 
annihilating with energies below 5 keV will be discussed later. It can be seen 
that the peaks of these spectra lie in the region 1/2 < rj < 1. This effect is due 
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to a pile-up of those gamma- rays from the annihilation of relativistic positrons 
which are emitted in the backward direction in the c.m.s. 

In Fig. 8 , this flux is compared with gamma-ray spectra from neutral 
pi-meson production from galactic cosmic- ray collisions (Stecker, 1967) and 
with the positron's own bremsstrahlung spectrum ( I n (E y ) and I B (E y ) re- 
spectively). It should be noted that these spectra are rigidly related since the 
ultimate source of I n (E y ) is the same cosmic-ray collision process generating 
the positrons and I B (E y ) is completely determined by I + (y) and n e . The total 
galactic bremsstrahlung spectrum is of course determined by the sum of the 
cosmic-ray positron and electron fluxes in the galaxy and may be expected to be 
at least twice as large as the I B ( 77 ) flux shown in Fig. 8 . We have also shown 
the expected bremsstrahlung gamma-ray flux from the observed cosmic-ray 
electron spectrum. 

Recent measurements have indicated that for 77 <2 the observed isotropic 
gamma- ray spectrum follows a power law of the form I obs (77) ~ if 2 3 (See 
review paper by Gould, 1967). An extrapolation of the observed isotropic 
spectrum is also plotted in Fig. 8 along with the conjectured extrapolation of 
Shen and Berkey (1968). Various authors have suggested that this flux may be 
due to Compton collisions between intergalactic positrons and electrons and the 
universal thermal microwave radiation (Hoyle, 1965; Gould, 1965; Felten, 1965; 
Fazio, Stecker and Wright, 1966; Felten and Morrison, 1966 or background 
from external galaxies (Gould and Burbidge, 1963; Silk, 1968). If this is indeed 
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the case, it may place severe restrictions on observations of the galactic 
AGS. 

IV. THE ANNIHILATION-GAMMA-RAY SPECTRUM FROM COSMIC-RAY 

POSITRONS ANNIHILATING NEAR REST 

We now come to the problem of determining the AGS from cosmic-ray 
positrons annihilating near rest. Because some aspects of this determination 
seem deceptively simple, there has been a tendency to oversimplify this problem 
in the literature. For this reason, I will first list various aspects of the prob- 
lem essential to an accurate analysis before proceding to treat them. 

A. In considering annihilations near rest, one must consider the pos- 
sibility of the intermediate formation of the bound electron-positron system, 
i.e., the positronium atom. At low energies, the cross-section for positronium 
formation becomes much greater than the cross-section for free annihilation. 

B. Once formed in interstellar space, a positronium atom will annihilate 
75% of the time into three photons. This situation contrasts sharply with the 
case of free annihilations where three-photon annihilations occur with a proba- 
bility of less than 1/2%. Therefore, the three-photon annihilation process, which 
produces a continuous spectrum from 0-0.5 MeV, must be considered along with 
the two- photon line- annihilations. Their relative importance depends directly on 
the fraction of positrons which ultimately form positronium. 

C. The most important source of the cosmic-ray positrons having energies 
greater than a few MeV is the source we have been considering, the decay of 
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secondary charged pi-mesons. There are, however, sources of relatively low 
energy positrons (less than a few MeV) which, as we shall see, may have a 
greater probability of being trapped inside the galaxy until they annihilate near 
rest. These sources are the beta- emitters, which are produced predominantly 
by low energy cosmic- ray interactions involving carbon, nitrogen and oxygen. 
Therefore, any observable 0.5 MeV line radiation from the galaxy may be pri- 
marily due to these beta emitters and that an observation of the intensity of the 
0.5 MeV line could supply information on the intensity of low-energy cosmic- 
radiation in the galaxy. 


V. GALACTIC POSITRONIUM FORMATION 

The cross-section for positronium formation by fast positrons in atomic 
hydrogen has been calculated by Cheshire (1964) and will be used as an approxi- 
mation for the interstellar medium. The ratio of positronium formation to free 
annihilation is only significant at non- relativistic energies and may be approx- 
imated by 


r 


p o s 


(r) 


CT A O') 


S(7) 


0 


10"6 ly 3 

10“ 2 (7 - 1)“ 2 
6.4 x 10 11 (7- 1) 
0 


for (7"1) > 10“ 2 

10" 4 < (7- 1) < 10“ 2 
2.5 x 10“ 5 < (7- 1) < 10“ 4 (34) 
10“ ® < (7-1) <2.5 x 10“ 5 
(7~1) < 10“ 5 


At energies of the order of the hydrogen binding energy, the probability 
-for positronium formation may be determined as follows (Deutsch, 1953): 


17 



Whereas free electrons and positrons may readily combine to form positronium 
with a binding energy just half that of hydrogen (the effective Bohr radius being 
2a 0 ), a positron needs a certain threshold energy in order to break apart an atom 
so that it may combine with its electron to form positronium. The binding energy 
of hydrogen is 13.6 eV; the binding energy of positronium is 6.8 eV. If we con- 
sider a fast positron impinging on a hydrogen medium and being slowed down by 
ionizing collisions, and if we assume that the last ionizing collision is equally 
likely to leave a positron with any kinetic energy between zero and Vj (V i being 
the ionization energy of the medium), then the fraction 

W i ~ 6.8 

f u * V. (35) 

1 

of positrons may form positronium. 

We must also take into account the possibility of inelastic collisions rapidly 
reducing the energy of the positron, not through ionization of the atom but through 
excitation. Therefore, if the first excitation level of the atom is V 1( and if 
Vj > (Vj - 6.8), then the fraction 

Y 1 + 6.8 - V. 

U * V L < 36 ) 

i 

should be a lower limit for the probability of positron formation. For a hydrogen 
medium 

V. = 13.6 eV and V x = 10.2eV, (37) 
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and therefore the probability of positronium formation, Pj], is 


0.25 < Pq < 0.50 . (38) 

For practical purposes, we may therefore use 

Pn = 1 ' (39) 

By a similar argument, we can show that the probability of the formation of 
positronium in its first excited state (n = 2, binding energy of only 1.7 eV) is 
quite small. 

VI. GAMMA-RAY SPECTRA FROM TWO- AND THREE- PHOTON 
ANNIHILATIONS AT REST 

The cross-section for free annihilation into £ photons, i.e., for the 
process 




(40) 


is of the order 


°A t £7 




a(^ 2 > 


°"A,2y ’ 


(41) 
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where 


a 



(42) 


Since we have seen that the positronium formation process should also be 
considered, we must also consider reactions of the type 

e + + e" - n (43) 

where the capital pi symbol will stand for positronium. Here again, we can 
neglect processes involving £ >3. 

The positronium annihilations, unlike the free annihilations, obey selection 
rules since positronium is in an eigenstate of charge conjugation C, and C is 
conserved in electromagnetic interactions. 

If £ photons are produced in the final state, then 

C = (-l) c . (44) 

It can be shown that under particle exchange, positronium obeys a kind of 
generalized Pauli principle and changes sign. The exchange of particles involves 


* Positrons may annihilate with bound electrons producing a single photon, with momentum being 
conserved by the binding nucleus. However, the cross section for this process is always less 
than zV(137)^, where the atomic number Z is almost always 1 or 2 under astrophysical con- 

ditions. Therefore, the single-photon annihilation can be neglected in the following discussion 


(Heitler, I960). 
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the exchange of the product of the space, spin, and charge- conjugation parts of 
the wave function. Therefore, 

(-l/ (-1) S + 1 C = -1 , (45) 

where •£ is the orbital angular- momentum quantum number and S is the spin 
quantum number. By combining Eqs. (44) and (45), we obtain the selection rule 

l + S = £ . (46) 


We may therefore specify the processes given by Eq. (43) as 


e 


+ 


+ e 


n ('s 0 ) 

L 2y . 


and 


(47) 


e + + e- - npSj) (48) 

^ 3y . 

Processes involving -£ f 0 can be neglected (Deutsch, 1953). 

The lifetime for positronium annihilation into two gamma-rays is given by 

T 2y = 4 ( 7 A, 2y v l'/'(°)l 2 ’ 



r 2 c 
r o e 

2a 0 3 n 3 


a 6 c 

2r 0 n 3 


(49) 
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or 


t. - 1.25 x 10~ 10 n 3 sec , 

V 2 ^ 

where n is the principal quantum number of the positronium state, a Q is the Bohr 
radius given by 


a o - a 2 r o ’ ( 50 ) 

v is the electron velocity, and |<//(0)| is the absolute magnitude of the wave 
function at the origin of the system. The effective Bohr radius for positronium 
is twice that of hydrogen, and this is taken into account. The factor of 4 in 
Eq. (49) takes into account the fact that for the two-photon decay of positronium, 
we know we are concerned with the "singlet" state, whereas in the case of free 
annihilation the probability of a singlet interaction (spins antiparallel) is 1/4 
and the probability of a triplet interaction (spins parallel) is 3/4. The factor 
1 1 //( 0 ) | 2 is just the effective electron density seen by the positron. 

The lifetime for positronium annihilation into three gamma rays is given 
by (Ore and Powell, 1949) 


or 



3y 


9v 


(-rr 2 - 9) 


7 

a' c 


r o n ' 


(51) 


t - 1.4 x 10 7 n 3 sec . (52) 

d S x -3y 
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The lifetimes for both the two- and three-photon annihilations of positronium 
are therefore so short that for considerations of galactic gamma-ray production, 
we may consider the annihilation to take place effectively at the time that posi- 
tronium is formed. Therefore, we conclude that 3/4 of the positronium formed in 
the galaxy annihilates into three-photons. 

The energy spectrum of the two-photon annihilation in the center-of- 
momentum system (c.m.s.) of the electron-positron pair is single line at 
E y - 0.51 MeV (rj = 1), as can easily be seen from considerations of conserva- 
tion of energy and momentum. The natural width of this line due to the uncertainty 
principle is small, being on the order of 

AE = ^ “ 5.3 x 10~ 12 MeV . (53) 

Dominant broadening can be expected to be due to the Doppler effect, and 
is of the order 


AE 

E 


AE 

Me 2 


■ 


(54) 


The effect of astrophysical conditions on the broadening of the 0.51- MeV 
line from two-photon annihilation can be determined as follows: 

In free e + - e" annihilations, we may consider a gas or plasma at tem- 
perature T. Then the distribution of the component of particle velocity along 
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the line of sight of the observed gamma ray is of the form 


f (/d|,) d/3 h 



-b/3 

e 


2 

II 


d/3,| 


( 55 ) 


where 


= Me 2 

b " 2kT ' (56) 

It then follows from Eqs. (23) and (36) that the spectral -line shape has the 
Gaussian form 


f(p) dr] 



[(T 1 


dr] 


(57) 


Equation (57) has the same form as Eq. (55), since the number of collisions in- 
volving velocity v is proportional to <xv -cr 0 c independent of v. 

The order of magnitude of the broadening is of the order of 


= (t) - »-‘ /2 - 1.8 x 10"'T« , (5; 

so that for T = 100°K, AE - 10 _1 keV. The Doppler width for two-photon posi- 
tronium annihilation will be discussed later. 

The energy spectrum from the three-photon annihilation is continuous, 
as allowed by conservation of momentum. It has been calculated by Ore and 
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Powell (1949) to be of the form 


. _ 2 fVl-T,) 2(l-77) 2 , ,, n 


2 -77 , 2(1-77) 


7? 


7 ] 


*11( 1 - 77) , (59) 




The function F(7j) is shown in Fig. 5; F(? 7 ) is normalized so that 


F(t?) dr) - 1 . (60) 

Jo 

VH. THE NUMBER OF COSMIC-RAY POSITRONS ANNIHILATING NEAR REST 

We now turn our discussion to the problem of the number of cosmic-ray 
positrons annihilating near rest. These positrons most likely come from two 
source: 

A) Positrons from the decay of secondary charged pi-mesons which were 
created at low enough energies to be trapped for a sufficiently long time in the 
galaxy to lose essentially all their energy before either annihilating in flight or 
escaping from the galaxy. 

B) Positrons from the decay of beta-emitting nuclei formed in collisions 
of low-energy cosmic-rays involving nuclei of carbon, nitrogen and oxygen. 
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The fraction of the original positron flux from the decay of secondary 
pi-mesons which annihilate near rest, f + , is given by 


f + 


r 

r dw 

I dz q + ( z ) exp 

Ji lr(w)|r s (w) 

Jl 


dz q + (z) 



( 61 ) 


This fraction was calculated numerically using Eqs. (24) - (29) and is given in 
Table 2 and Fig. 6 for various possible mean leakage times, T^. Table 3 
gives the corresponding values of Q rest the total number of positrons from 
pi-meson decay per cm 3 per second annihilating below 5 keV for the halo model 
of Table 1 (Q rest „• being defined by the numerator of Eq. (61) and the resulting 
gamma-ray fluxes). 

In the case of an infinite leakage time (all positrons being trapped and an- 
nihilating in the galaxy) we find that 80% of the positrons produced annihilate 
near rest, a figure which is in perfect agreement with that given by Heitler 
(1960) as an asymptotic value for the annihilation of ultrarelativistic positrons 
when the dominant energy loss comes from ionization. However, for the 
leakage time usually considered as plausible for the galactic halo, 10 8 years 
(corresponding to x - 5 gm/cm 2 ), only 1-2% of these positrons annihilate near 
rest. 
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We next consider the positrons produced in the decay of beta-emitting 
nuclei produced in low-energy cosmic -ray collisions. The important reactions 
to be considered are listed in Table 3 (Audouze, et al., 1967). Whereas the 
majority of pi-mesons are produced by cosmic-rays having energies over 500 
MeV, beta -emitting nuclei can be produced by cosmic-rays with energies down 
to below 20 MeV. The positrons produced have energies less than 5 MeV with 
the singular exception of those produced from the decay of N 12 . Using the cross- 
section data given in Table 3, and the quiet-sun observations of the cosmic-ray 
spectrum from 20-1000 MeV/nucleon (Comstock, et. al., 1966). We obtain the 
value for the fluxes given in Table 4. This value probably represents a lower 
limit to the true astrophysical value since over 95% of these beta-ray positrons 
will remain trapped in the galaxy and annihilate near rest. For mean leakage 
times less than about 400 million years most of the positrons annihilating near 
rest come from the beta-emitting unstable nuclei produced in cosmic-ray p-C,N,0 
interactions. 

Equation (34) may now be used to determine the amount of positronium 
being formed by positrons combining with electrons as a function of energy. 

The result of this calculation is shown in Fig. 7. This figure shows the percentage 
of positrons which, after having survived to reach a 10 keV kinetic energy, sur- 
vive to reach lower energies. The dashed line indicates the survival fraction 
for free annihilation only; the solid line indicates the survival fraction when 
positronium formation is taken into account. Figure 5 shows that almost all of 
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the positrons annihilating near rest do so through intermediate positronium 
formation with an average energy of about 35 keV. Therefore 75% of the 
resultant annihilations occur via the three-photon channel and 25% produce a 
two-photon line with a Doppler width of about 5 keV. 
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Table 1 


Average Astrophysical Parameters for the Galactic Disc and Halo. 



n e (cm -3 ) 

HOug) 

P y (e V/cm 3 ) 

T£ (10 6 yrs.) 

L eff (cm) 

Disc 

0.5-1 

6 

0.45 

2-4 

3 x 10 22 


lira 

CO 

0.65 

100-300 

5 x 10 22 


Table 2 

Fraction of Positrons from 7T + -decay which Survive to 
Annihilate Near Rest (T + < 5 keV). 


T^ (10 6 yrs.) 

U (%> 

00 

80 

600 

20 

300 

9.3 

100 

1.6 

60 

0.61 

30 

0.14 

10 

0.022 

6 

0.0065 

3 

0.0025 
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Table 3 


Principal Reactions Leading to Production of /3 - Emitting Nuclei 
(See Audouze, et al. (1967)). 


Reaction 

a (mb) 

— 

Decay 

Halflife 

Positron Energy (MeV) 

C 1 2 (p,3p2n)B 8 

? 

B 8 (/3 + )Be 8 

0.78 sec 

1.4 

C 1 2 (p,p2n)C 10 

-3 




N 1 4 (p,2p3n)C i0 

? 

C 10 (/3 + )B 10 

19 sec 

1.9 

0 16 (p,3p,4n)C 10 

<10 




C 1 2 (p,pn)C 1 1 

-70 




N 1 4 (p,2p2n)C n 

-30 

C 1 1 (/3 + )B 11 

20.5 min 

0.97 

0 1 6 (p,3p3n)C n 

-10 




C 1 2 (p,n)N 1 2 

? 




N 14 (p,p2n)N 12 

? 

N 12 (/3 + )C 1 2 

0.011 sec 

16.4 

0 1 6 (p,2p3n)N 12 

? 




N 1 4 (p,pn)N 1 3 

-15 

N 13 (/3 + )C 1 3 

10.0 min 

1.19 

0 16 (p,2p2n)N 13 

-10 




N 14 (p,n)0 14 

-50* 

0 14 (/3 + )N 14 

71 sec 

1.18,4.14 

O 1 6 (p,p2n)0 14 

<10 




0 16 (p,pn)0 15 

-50 

0 15 (/3 + )N 15 

2.06 min 

1.74 


? - Unknown and not included in estimate of positron production. 

*This cross section has a value of about 100 mb at energies below 12 MeV but 
is negligible above 150 MeV. 
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Table 4 


Positron Annihilation Rate and Gamma-Ray Flux from 
Positrons Annihilating Near Rest. 


Source 

Q-T.rest Sec_1 ) 

Flux* (2-y-line) 
(cm -2 sec -1 sr -1 ) 

Flux* (3y-continuum) 
(cm -2 sec -1 sr _1 ) 

7T + 

T \ (10 6 yrs.) 





00 

1.4 X 10" 27 

2 x 10" 5 

1 x 10" 4 


600 

3.5 x 10“ 28 

5 x 10" 6 

2.5 x 10 " s 


300 

1.6 x 10~ 28 

2.5 x 10" 6 

1 x 10 -5 


100 

3 x 10" 29 

5 x 10" 7 

2 x 10“ 6 


60 

1 x 10' 29 

1.5 x 10" 7 

7 x 10" 7 


30 

2.5 x 10" 3 ° 

4 x 10" 8 

2 x 10" 7 


10 

3.8 x 10" 31 

6 x 10" 9 

2.5 x 10“ 8 


6 

1.1 x 10" 31 

1.5 x 10" 9 

i 

7.5 x 10" 9 

’ 

3 

4.4 x 10‘ 32 

7 x 10' 10 

3 x 10" 9 

p-C,N,0 

2.5 x 10“ 28 

4 x 10" 6 

1.6 x 10" s 


* Assuming (nL eff ^ = 3 x 10 21 cm -2 . 
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min 



0 1 2 3 4 5 


E y /m e c 2 (17) 

Figure 1. The kinematic limits on the positron Lorentz factor involved in the determination 
of the laboratory-energies of the annihilation-gamma-rays produced. 
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Figure 2. The source spectrum of positrons produced from the decay of the positive pi-mesons 
formed in cosmic-ray collisions (Ramaty and Lingenfelter; 1967). 
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T ^ = 00 \ yX=30 g/cm 2 T^=6xl0 8 yr 



0 12 3 4 


log y 

Figure 3. Various positron equilibrium fluxes for the halo model of the galaxy given for 
various approximate mean path lengths and mean leakage times. 
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Figure 5. 
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V 


Energy spectrum of gamma rays resulting from the three-photon annihilation of an 
electron and a positron (from Ore and Powell, 1949). 
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SURVIVAL PERCENTAGE 



Figure 7. The percentage of positrons which, after having survived to reach an energy of 
10 keV, survive to reach lower energies. The dashed line indicates the survival fraction 
found by taking into account free annihilation only; the solid line indicates the survival 
fraction when positronium formation is also taken into account. 
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Figure 8. Isotropic gamma-ray spectra from various astrophysical processes. 
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